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Abstract: In this talk, I will argue that the HERA experimental data show that the 
typical parameter (k) responsible for the value of the shadowing corrections (SC) in DIS is 
so large that the BFKL Pomeron is hidden under SC. The SC turn out to be large enough 
but mostly for the gluon structure function which is not well determined by the available 
experimental data and by the current theoretical procedure. 



In this talk I am going to answer two questions: 

Ql: Where is the BFKL [| Pomeron? 

Q2: Where are shadowing corrections (SC)? 

Actually, the answers have been presented in our paper |]], but here I will discuss them in 
more details. 

First, let me explain why it is reasonable to ask such questions. Indeed, at first sight, the 
situation looks very transparent, namely, the HERA data can be described by means of 
the usual DGLAP evolution equations without any other ingredients such as the BFKL 
Pomeron and / or SC ( see any of plenary talks during the past three years). My personal 
opinion is that this fact brought more questions than answers since we need to show ( to 
justify theoretically our approach) that the corrections due to the BFKL dynamics and/or 
due to the SC are negligible small at least at the HERA kinematic region. If it is not so 
( as I will show below) the DGLAP approach is not better or worse than any other model 
developed to describe the experimental data. The main goal of this talk is to show that the 
experimental data from HERA confirm that both the BFKL contribution and the SC should 
be rather large in the HERA kinematic region. 

Actually, everything that I want to tell is given in Fig.l, but I need to explain what are 
plotted in this figure. 

1. < 7 > = \ and < 7 > = 1. 

Let me recall a standard procedure of solving of the DGLAP evolution equations. 
The first step: we introduce moments of the structure function, namely, 

xG(x,Q 2 ) = — [ e~ u Hl/x) M(lu, Q 2 ) du, 

2,711 JC 

where contour C is located to the right of all singularities of moment M(u, Q 2 ). 
The second step: we find the solution to the DGLAP equation for moment 

dM(uj,Q 2 ) 



(ilnQ 2 

The solution is 



7 (cu)M(u;,Q i 



M{lu,Q 2 ) = M(u,Ql) ■ e ^ )ln(Q /Q ° } . (2) 

Here M(cu, Qq) is the nonperturbative input which should be taken from experimental data 
or from "soft" phenomenology ( model). 

The third step: we find the solution for the parton structure function using the inverse 
transform, namely: 

xG(x,Q 2 ) = f — e -MiA) + 7Hin(Q 2 /Qg) M (uj,Q 2 ) . (3) 
Jc 2rri 
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Therefore, to find a solution of the DGLAP equation we need to know the nonperturbative 
input M(u,Qq) and the anomalous dimension 7(0;), which we can calculate in perturba- 
tive QCD. The anomalous dimension 7(0;) has been calculated in pQCD and the result of 
calculations can be written in the form: 

/ \ BFKL/ a sN c a s N c 2 

7(w) = 7 ( ) + <*s7iM + a s l2(uj) , (4) 

7r u iru 

where both functions 71 and 72 are known as well as ^ BFKL . Using Eq. (|j) we can discuss 
what has been done in the global fits [B]. The value of the anomalous dimension has been 
calculated in as and a 2 s orders ( two last terms in Eq.(5)) and the nonperturbative input 
has been taken in the form M(u, Qq) oc with ujq ~ 0.2 - 0.3. This means that the 

structure function at Q 2 = Ql increases as x~^° at x — > f\. However, one can see that the 
ryBFKL ghouid Iqq essential in the region of low x where u — > since 



^BFKL, a sN c a s N c _ a s N c 1 u - u L 

7TC<J TIL) ^4 7TC<J 2 V ^ 

where cj^ and A have been calculated 0. 

This equation reflects the main properties of the BFKL Pomeron: the limited value of 
the anomalous dimension and the importance of all terms of the order of (2fi^s) n in the 
region of small u. All attempts to estimate the values of the BFKL terms in the anomalous 
dimension || show that they are essential in the HERA kinematic region. Here, we choose a 
different way of presentation of this well known fact, namely, we introduce average anomalous 
dimension < 7 > which is equal to 

1 dxG(Q\x) Jc jfeg e" ln ( 1 /") + ^) l ^ 2 /^M(co,Q 2 ) 

< 7 > ~ xG{Q 2 ,x) ' d\n{Q 2 /Ql) ~ J c ^ & HV*)+rt«)WQ a /<%) M(w,Ql) ' ( } 

Function < 7 > describes the behaviour of the anomalous dimension quite well since at low 
x the deep inelastic structure functions can be calculated in the semiclassical approach 
in which, for example xG(Q 2 , x), is equal to 

xG(Q>,x) = C(ln(l/x),ln(QVQS)) ■ x< ^J Mfl , m)> ■ (g)«*^>*«W«»» (7) 

where functions C,< u > and < 7 > are smooth function of ln(l/x) and \n(Q 2 / Ql) . 

In Fig.l we plotted two lines with < 7 > = | and < 7 > = 1 . We expect a large the 
BFKL contribution in the kinematic region between these two lines and one can see in Fig.l 
that we have penetrated this region at HERA. 



* Strictly speaking this statement is correct for two global fits: MRS and CTEQ. The GRV fit has a 
different initial condition, namely, the evolution has been started at very low value of Q 2 but with the initial 
distribution which is flat at low x . 
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Figure 1: Contours for < 7 > = 1 and 1/2 and k = 0.6,1,1.6 for the GRV95 gluon density 
and HERA kinematic region. 
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2. K. 

From HERA data we can evaluate also the probability k of the parton - parton (gluon - gluon) 
interaction, which is given by M, M 



2 .a(GG) 3 nag „ 2 \ 



where xG(x, Q 2 ) is the number of partons ( gluons) in the parton cascade and R 2 is the 
radius of the area populated by gluons in a nucleon. a{GG) is the gluon cross section inside 
the parton cascade and was evaluated in @. 

The observation is that we know from the HERA data both the value of the gluon 
structure function and the value of R 2 in Eq. (|8]). Indeed, the available parameterizations 
such as MRS, CTEQ and GRV f| give the value of the gluon structure function with 
sufficiently large differences in its value. However, this difference is less than 50% and 
becomes smaller with improvement of the experimental data on F2(x,Q 2 ) ( see Fig. 2 ). 

The most important and new information is the fact that using HERA data on photo- 
production of J/\l/ meson the value of R 2 can be estimated as R 2 < 5GeV~ 2 0. 

Indeed, (i) the experimental values for the slopes ( see Fig.3 ) are B e \ = AGeV~ 2 and 
Bi n = 1-66 GeV~ 2 and (ii) the cross section for J/\l/ production with and without proton 
dissociation are equal ||. Taking into account both facts we can estimate the value of R 2 ( 
see Ref.[[J for details) which appears in calculation of the SC (Glauber corrections) due to 
integration over the momentum transferred (q]_ = \t\) along the gluon ladders (see Fig. 4) 
neglecting t dependence of the upper vertex in Fig.4. 

In Fig.3 we show the picture for the diffractive production of J/^f in the additive quark 
model, in which two radii naturally appear as the radius of the hadron and a proper radius 
of the constituent quark. In all our estimates we did not need this particular model but it 
is interesting to mention that our estimates give the same value of the average radius as in 
the additive quark model. 

It should be stressed that such an estimate gives the value for R 2 which lead to the value 
of the cross section for the double parton scattering measured by the CDF collaboration at 



the Tevatron 10 



Let us discuss this point a little bit in more details. The CDF collaboration measured 
the processes of unclusive production of two pairs of "hard" jets with almost compensated 
transverse momenta in each pair and with almost the same values of rapidities. Such pairs 
can be produced only due to double parton collision and their cross section can be calculated 
using the Mueller diagram given in Fig. 5. 

The value of the double parton scattering cross section can be written in the form ( see 
Fig.5)|g: 

a DP = m , (9) 

2(J eff 
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Figure 2: The value of the gluon structure function in different parameterizations 



5 




6 



Figure 5: Inclusive production of two pair of "hard" jets in the double parton scattering. 



where, for simplicity, we consider the production of two ( QiQi and Q2Q2 ) quark - antiquark 
pairs. Factor m in Eq. (Q) is equal to 2 for different quarks ( Q% 7^ Q2) and to 1 for identical 



quarks. The value for a e ff is measured to be 14.5 ± 1.7 ± 2.3 mb. Our estimates [|TTJ for 
diagrams of Fig.5 using the two radii picture give cr e // ~ 17mb. It means that the effective 
radius R 2 = 5GeV~ 2 could be even overestimated. 

Using the GRV parameterization for the gluon structure function and the value of R 2 = 
5GeV~ 2 , we obtain that k reaches 1 at HERA kinematic region ( see Fig.l ), meaning 
shadowing corrections should not be neglected. In Fig.l we plotted three curves with values 
of k equal to 1.6, 1 and 0.6, respectively, to illustrate a possible range of k using CTEQ and 
MRS parameterizations. 



Al: 

The answer to the first question one can read from Fig.l. Indeed, the kinematic region where 
the BFKL Pomeron ( the BFKL corrections to the anomalous dimension ) could be sizeable, 
namely, the region between curves < 7 > = 1/2 and 7 = 1 , is located to the left of the curve 
with k = 1 where the SC should be essential. Therefore, we can conclude that the BFKL 
Pomeron is hidden under large SC and cannot be observed. To illustrate this point and to 
show what is the influence of the SC on the behaviour of the average anomalous dimension 
we plotted in Fig.6 < 7 > given by Eq.(6) but using the Glauber - Mueller formula for the 
SC for the gluon structure function ( see Ref. |1| ). The calculations were performed for the 
gluon structure function at fixed impact parameter (bt = 0), where we take 

xg(x,Q 2 ) = J d 2 b t S(b t )xG(x,Q 2 ) , 

with 

1 "l 
S(b t ) = —e~W . 
ttR z 
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One can see, that the average anomalous dimension turns out to be smaller that < 7 > = ~. 
Therefore, the BFKL Pomeron will not be seen even if we will take the SC at the minimal 
rate given the Glauber - Mueller formula. 

A2: 

It turns out that the SC is not very big for F 2 (Q 2 ,x) which has been measured experimen- 
tally ( see Ref. Jl| for details). However, the SC for xG(Q 2 ,x) should be large. To illustrate 
this point we plot in Fig. 6 the ratio Ri = xG(q 2 ,x) Grv calculated i n Ref- flTJ. Comparing 
Fig. 6 with the value of xG(Q 2 , x) in current parameterizations we can conclude that in spite 
of sufficiently large SC our knowledge of the value of the gluon structure function is so poor 
that we can absorb all SC in the uncertainties of its value. For example, taking into account 
the SC the GRV gluon structure function will be able to describe the new experimental data 
while without the SC the GRV parameterization has been ruled out by experiment ( see the 
last picture in Fig.2.). 




y=ln(1/x) 

Figure 6: Average anomalous dimension, given in Eq.(6) for the SC taken in Mueller- 
Glauber formula for b t = and ratio R\ for the gluon structure function. 

It should be stressed, that our answers bring several problems that have to be solved in 
the nearest future. These are three of them: 

PI: 

It is very likely that the BFKL Pomeron will be hidden under SC not only for the deep 
inelastic proton structure function but in all specially invented processes to extract the BFKL 
Pomeron since the typical size ( the value of R 2 ) is smaller in all such processes than in the 
case of the deep inelastic structure function. 

P2: 
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The fact that the SC will take place of the BFKL Pomeron does not mean that the 
the experimental cross section will be the same as in the DGLAP evolution equations or 
in the Monte Carlo simulations based on the DGLAP evolution. The difference should be 
calculated to be discussed. 

P3: 

The size of the SC for the gluon structure function crucially depends on the initial gluon 
distribution. In our estimates we pretended that the GRV parameterization guessed correctly 
this initial distribution. The only argument is the fact that the GRV parameterization 
describes the experimental F 2 (x,Q 2 ) at small Q 2 ~ lGeV 2 . 

Alternative answers: 

AA1: 

The first alternative answer has been proposed by R. Thorne (see Ref. ]12|] ) who demon- 
strated that the correct inclusion of the BFKL anomalous dimension allows to improve the 
comparison with the experimental data. However, Fig. 7 shows that the value of the gluon 
structure function extracted from the experimental data is not very different from the pre- 
vious analysis without the BFKL anomalous dimension. Therefore, the experimental data, 
perhaps, does not contradict the existence of the BFKL Pomeron, but this fact cannot change 
considerably the value of k. It means, that the value of the SC is still big even in the analysis 
taking into account the BFKL Pomeron ( see KMS paper for details). In Fig. 7 is given 
some next to leading order corrections to the BFKL Pomeron. One can see that they are 
essential and diminish the value of the gluon structure function but still not more than in 
two times, which we evaluate as a typical error in the value of extracted gluon structure 
function. 

AA2: 

The BFKL Pomeron is not seen in the data because the next to leading correction is 
essential and they change crucially the main properties of the BFKL Pomeron. Fortunately, 



the next order correction to the BFKL Pomeron (NOBFKL) has been calculated [13| and 



the community of experts has started to understand the influence of the NOBFKL [14 



[ p~5|j on the value of the gluon structure function. It turns out the the NOBFKL Pomeron 
has a much smaller intercept or in other word the power - like behaviour, namely, xG * 
x,Q 2 ) oc (—) UL , still remains but ujf° = uj^°(1 — 3.4 «s). It means that the energy ( 
x) dependence becomes milder and it makes the NOBFKL Pomeron not so pronounced as 
it was in the leading order. Nevertheless, the first numerical estimates show that the value 
of the gluon structure function changes but not significally. Indeed, Fig. 8 that was taken 



from Ref. |15| shows that the gluon structure function in the analysis with the NOBFKL 
anomalous dimension typically on 30% less than the gluon structure function without the 
BFKL contribution. 

My conclusions: 
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Figure 7: The gluon structure function with and without the BFKL Pomeron ( picture is 



taken from the KMS paper ||12|| ). 
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Figure 8: The gluon structure function with and without the NOBFKL Pomeron ( picture 
is taken from the Blumlein and Vogt paper 
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These two examples of the alternative answer show that (i) it is not so easy to change 
significally the value of the gluon structure function and diminish it more than in two times; 
and (ii) the real accuracy of the value of the gluon structure function is rather big, about 
50%, in spite of the fact that the difference between two sets of gluon structure functions 
( so called global fits: MRS and CTEQ ) became much smaller using new more accurate 
experimental data (see Fig.2). Our errors are mostly theoretical ones. In my opinion, we 
cannot change ( diminish ) the value of parameter k and therefore, accordingly to Fig.l, we 
have to deal first with the SC and only after that to take into account the BFKL Pomeron 
with all possible corrections. 

Finally, I would like to emphasize that you got my personal answers. Perhaps, you 
have different ones. The only point, which I insist on, is that, before answering these two 
questions, we cannot trust the DGLAP evolution more than any other model. 
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